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A recent argument of Oerding shows that our calculation of the quantity A, which determines 
the amplitude of the asymptotic decay of the particle density in 2 < d < 4, was in error. Instead it 
is simply given by A = no, the initial density, for uncorrelated initial conditions. 
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In recent years much progress has been made in ap- 
plying field-theoretic renormalization group methods to 
models of diffusion- limited chemical reactions. One such 
example was a study by the present authors of the 
A -\- B ^ ^ annihilation reaction with random initial 
conditions and equal particle densities Our purpose 
here is to correct a mistake contained in this work, and 
to revise the conclusions accordingly. In particular, we 
now find for spatial dimensions 2 < d < 4 that the A, B 
particle densities decay in time with the universal form 



{ait)) = m) 
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where (a(0)) = {b{0)) — uq, and D is the diffusion con- 
stant. This amplitude is exactly that obtained in the ap- 
proximate calculation of Toussaint and Wilczek |^ , and 
the initial density dependence is the same as found by 
Bramson and Lebowitz in the limit of an instantaneous 
reaction In our previous study we identified nonuni- 

versal corrections to the amplitude in (|l]) of order ng^^ 
[1|. However, in analogy with an argument of Oerding 
[i| for the same reaction in a random shear flow, it can 
be shown quite generally that the sum of all such cor- 
rections in fact vanishes. We now briefly paraphrase his 
arguments in the context of the present problem. 

In general one begins with a master equation descrip- 
tion of the diffusive and reactive particle dynamics. This 
is mapped to a field theory using by now standard tech- 
niques |l|,p|,D- It is convenient to replace the fields a(x, t) 
and b{x,t) with their sum (0) and difference {ip), result- 
ing in the action H 



d'^xdt^<P{df 

+ A2( 
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(2) 



The relation between the couplings in (g) and the origi- 
nal master equation parameters is given in ref. 1. Note 
that the diffusion constant is absorbed into a rescaling of 



time, and the conjugate fields (j) and ip are introduced in 
the mapping. 

The renormalization group (RG) analysis reveals that 
the critical dimension for the coupling A2 is dc — 2, i.e. 
A2 is irrelevant for d > 2. Since the dynamic RG re- 
lates evolution in time to renormalization group fiows, 
one finds the asymptotic behavior of the theory is given 
for d > 2 by an effective action of the form 

S,s = J d''xdtS^^{dt-\/^)<l> + ij{dt-V^)i^ (3) 

Under the renormalization fiows for d > 2 one has 
Ai — > Aoff, and, more importantly, new initial terms are 
generated such as ^A-ip^. While many such terms may 
appear and be relevant for < 4, it was shown in ref. 
1 that for 2 < d < 4 the asymptotic behavior of the 
density and correlation functions is determined solely by 
A. Hence we now focus on the generation of this cou- 
pling. (For d > 4 one obtains the rate equation result 
(a(t)) cx 1/Kst.) 

Denote the sum of all diagrams in an expansion of (|) 
which terminate at time t with two external ip lines by 
4>{t)'^^{t). This is illustrated schematically in fig. 4(a) of 
ref. 1. Since H(i) is damped for times t ^ l/n^Ai, the 
replacement 



i^{t)^ii{t) ~ i>{o)^s{t) / dt'nit' 



(4) 



is valid for calculating asymptotic quantities, and serves 
to define A. 

The sum of all such diagrams containing no loops 
is given below in fig. 1(a), and yields a contribution 
A(°) = no, the initial density. Here, following the no- 
tation of ref. 1, the dashed line represents the ip propa- 
gator, the heavy solid line is the (j) classical (tree-level) 
response function, which is the (j) propagator dressed by 



1 



the initial density, and the wavy hne is the classical den- 
sity. Three- and four-point vertices have coupling con- 
stants Ai and A2, respectively, with signs which can be 
determined from (||). In ref. 1 we assumed that all other 
diagrams in n(t), except those which dressed the A2 ver- 
tex, could be accounted for by taking Ai ~> Xcs in the 
tree level diagrams. This left A*^"^ unchanged, and the 
corrections were found to yield an expansion in powers 
of uq of the form A — jiq — C'tIq^'^ + . . ., with C nonuni- 
versal (see equation (1.5) in ref. 1). However, Oerding 
has shown that A can be calculated exactly from the full 
action (0), with the result B 



implies that the contribution from fig. 1(c) is 



A = no. 



(5) 



Since A itself is not renormalized fyl this must be the 
correct value for the effective action (^, in contradiction 
with the previous result. 
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Hence A = A^"), which leads to §) §. 

Using the correct and universal result for A in the anal- 
ysis of ref. 1 yields the asymptotic density (|l|). In ref. 1 we 
also derived the correlation functions and the density in 
the case of unequal diffusion constants Da ^ Db- Since 
these results are expressed in terms of A, they may be 
regarded as correct and universal once the substitution 
(^) is made. However, we stress that these results are 
universal only for truly random initial conditions. Cor- 
relations initially present are expected to modify A in a 
non-universal way. 

We thank K. Oerding for bringing this error to our 
attention. 



(a) (b) (c) 

FIG. 1. Diagrams which contribute to n(t), hence to A. 
The symbols are defined in the text and in ref. 1. 

Oerding's result follows from observing that all dia- 
grams in H(t) belong to one of the three groups in fig. 
1. The shaded area represents the same set of diagrams 
in fig. 1(b) as in fig. 1(c). These may be attached to 
the highlighted vertex (represented by a dot) by either 
two i) or two 4) lines: in either case, the sign of the A2 
vertex in fig. 1(b) is the opposite of the Ai vertex in fig. 
1(c), as can be confirmed by S (||). Therefore, writing 
the contribution from fig. 1(b) as 



A^") = A, 



dt fit) 
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